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1.0 INTRODUCTION 


In remote sensing applications there is a need for an atmospheric 
model to simulate the response of a satellite sensor for realistic, var- 
iable atmospheric conditions. In addition, there is a need for an atmo- 
spheric model to "correct" the sensor response for expected or anticipated 
atmospheric conditions which exist for a particular region of interest. 

The remote sensing of the Earth's surface and the interpretation of 
various surface elements is a complex problem for many reasons. First, 
there is the variability of the surface itself; the reflectance and emis- 
sivity of surface features are dependent upon factors which are either dif- 
ficult to model or predict. Thus, the intrinsic radiance of a species of 
vegetation will depend upon the growing season, soil moisture, time of day, 
wind speed, precipitation history, and perhaps other factors which are dif- 
ficult to analyze. Second, atmospheric effects can play a role in varying 
the radiance of the surface material at the surface and also in the radi- 
ance along the path which connects the sensor and the target element. 
Thus, one can have a variable atmosphere as a result of the spatially and 
temporally varying aerosol or the variable gaseous components such as water 
vapor, ozone, and carbon dioxide. In addition, there is the presence of 
clouds and, therefore, cloud shadows which must be accounted for in the 
analysis of remote sensor data. 

Remote sensing means different things to different groups. One may 
be interested in estimating the proportion of crops or materials in a given 
region. Another may be interested in the identification of various fea- 
tures based upon spectral properties. In early work in remote sensing, 
investigators attempted to create signatures of various materials based 
upon their actual spectral characteristics in a known area. Then the sen- 
sor was transported to another area and each pixel of data was then tested 
to determine how it compared spectrally with the known signatures. By des- 
ignating a specific area for investigation one could then classify the 
pixels according to their spectral classes. Another method can be used 
however; one which is based upon models of various materials. The models 
can be general so that they are characteristic of many of the variables 



which are used to represent actual data. Thus, a particular crop model 
might have the capability of predicting spectral reflectance in terms of 
the period during the growing season, the time of day, soil moisture, and 
wind speed. This is then an attempt to remove some of the variability 
associated with various material classes. 

A similar problem exists with the atmosphere. In the extension of 
spectral signatures from one area to another the atmosphere will be dif- 
ferent and, therefore, errors will be introduced into the processing of 
images for classification analysis. A deterministic atmospheric radiation 
model would certainly be desirable. Unfortunately, it is impossible to 
predict the optical state of the atmosphere throughout the time span and 
spatial region over which remote sensing takes place. In lieu of such 
models, we shall investigate the statistics of basic atmospheric optical 
properties so that we can determine a measure of the variability of the 
atmosphere for particular times and areas of concern. We base our sto- 
chastic model on actual data for optical thickness as measured by sensors 
throughout the world. 
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2.0 MODEL PARAMETERS 


In remote sensing technology there are many parameters which are 
used in the description of surface features and the atmosphere. Because we 
are considering the atmosphere in this investigation rather than the de- 
tails of the surface, we will consider the atmospheric parameters in some 
detail. The primary parameters on which the radiation at the sensor 
depends are those listed in Table 1. 


Table 1. List of Primary Parameters 

for Remote Sensing 

Target Bidirectional Reflectance 

A A 

p(q,q') 

Background Albedo 

P B 

Spectral Optical Thickness 

V 

Target Directional Emissivit.y 

e{a) 

Background Emissivit.y 

e B 

Single Scattering Albedo 

w o 

Scattering Phase Function 

p(cosx) 

Density of Each of i-th Gas Component 

n. 

Extraterrestrial Solar Irradiance 

i 

E 

Geometry 

A , fi 

0 


Many of these primary parameters are, in turn, related to more basic quan- 
tities such as those listed in Table 2. These quantities are used in 
models for the specification of the primary parameters and should be con- 
sidered for the detailed modeling of the inherent features of surface and 
atmospheric quantities. It should be noted that we have not included the 
functional dependence of each parameter in the tables. It will be assumed 
that the physical (radiometric) and optical quantities are spectral and de- 
pend upon the wavelength of the radiation. 
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Table 2. List of Basic Quantities for Primary Parameters 

Barometric Pressure 

P 

Aerosol Index of Refraction 

m 

Aerosol Size Distribution 

n(r) 

Aerosol Shape Factor 

s 

Particulate Number Density 

N 

Atmospheric Temperature 

T a 

Cloud Temperature 

T c 

Target Temperature 

T t 

Background Temperature 

T b 

Water Vapor Pressure 

e 

Spectroscopic Parameters 

A i 


We now describe briefly the radiometric quantities which are used in 
the deterministic and stochastic atmospheric radiation models. The surface 
radiance, that is, the upward-directed radiance at a surface element is 
given by the following equation: 


i (n) = / n .n'p(n,n')L. (n')dfl'. 


in 


( 1 ) 


where n $ is the unit normal surface vector, p(Q,fT ) is the bidirectional 
surface reflectance which is dependent upon the directional unit vector for 

A A 

the incoming radiation o,' and the vector for the outgoing radiation Q,. The 

A 

incoming radiance L in (n ) , in turn, depends upon the atmosphere and the 
background surface albedo. It is the task of agronomists, hydrologists, 
geologists, and others to present specific models of the bidirectional 
reflectance for their respective disciplines. Once this is done the atmo- 
spheric scientist can use his models of the atmosphere to generate the 
downwelling radiation field so that one can calculate the surface radiance 
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L (n). For many examples one can use the assumption of a diffuse surface 
in which case the bidirectional reflectance is constant with respect to 
direction and we have for the surface radiance 



E T0T’ 


( 2 ) 


where p is the surface albedo and Ejqj is the total (solar plus diffuse 
sky) downward irradiance on the target. 

The radiance at a sensor aperture can be expressed as 


L " L o T + L p’ 


where T is the transmittance along the path between the target and the 
sensor and Lp is the path radiance. These two quantities are dependent 
upon the scattering and absorptive properties of the atmosphere, the 
geometry, and in the case of path radiance, the surface albedo. Many 
models exist which can be used to calculate the transmittance and path 
radiance, some of which are summarized by Turner [1]. Because of the 
relatively simple closed-form expressions used in a modified two-stream 
model, we will use the model of Turner [2,3,4] in our analysis. 


Because of the ease with which we could obtain data and because 
there is a reasonably good understanding of atmospheric particulates, we 
will concentrate on the stochastic nature of the atmospheric aerosol in 
this investigation. In a previous study, Turner considered the repre- 
sentation of the spectral optical thickness of the atmosphere in terms of a 
random variable. It should be recalled that the aerosol optical thickness 
of Earth's atmosphere can be expressed as the following: 




(x.z)dz. 


(4) 


where k ^ (x,z) is the volume extinction coefficient at altitude z. Because 
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the particulates are so variable we can consider t 0 ,^(x) as the most var- 
iable part of the optical thickness, at least in the spectral regions where 
gaseous absorption is weak. Thus, the aerosol optical thickness can have 
any positive value, i.e., 

0 ^ t q ^(x) < °° ( 5 ) 

Another quantity of considerable importance in radi ati ve -transfer 
theory is the single-scattering albedo w Q (X). It is the ratio of the 
volume scattering coefficient to the volume extinction coefficient and is a 
measure of the amount of scattering. It has the following range: 



3.0 


STOCHASTIC MODEL 


In this section we consider the formalism for the statistical 
properties of the highly variable atmospheric components. 

3.1 Sensor Response 

A typical sensor has detectors the output of which is usually in 
volts. In any case, we can write the sensor output 



L(X) S k (A)dX, 


(7) 


where $ k (x) is the spectral response function of the sensor for the k-th 
channel. If the radiance L(x) is known then the integration can be per- 
formed and the response s k can be determined. 

The response s k and the radiance L(x) are considered to be stochas- 
tic quantities i.e., they have properties which are associated with random 
variables. As stated earlier, the randomness arises from two factors; 1) 
the unknown or statistical nature of the surface properties, and; 2) the 
statistical nature of the atmosphere. Thus, we consider the mean and 
covariance of the sensor response for channels k and k' i.e., 

R * = E w 

C kk‘ = E {( s k" R k^ s k'" R k^} * 
where E denotes the expectation. 

3.2 The Covariance Matrix 

For multispectral remote sensing applications one can use the 
maximum likelihood decision rule for classification. Let us designate the 
logarithm of the radiance by a vector X^, i.e.. 


(8) 

(9) 
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x k = In L k . 


( 10 ) 


One can then use the n-dimensional normal distribution function i.e.. 


f(x 


l ,x 2 


,...,x n ) = 


(2tt f /Z JT 


exp { 


k,k -1 


(ID 


when C is the determinant of the covariance matrix C kk ' * We will now 
develop an expression for the covariance matrix for the spectral radiance 
in terms of the variances, means, and covariances of the known quantities. 
Let us consider the stochastic parameters to be optical thickness x, 
single-scattering albedo w, and the background albedo p g . The spectral 
radiance for the i-th pixel and the k-th channel can be represented ap- 
proximately by 

L ik (T,u),p B ) = P ik F k (x,w,p B ) + G k (x ,w ,p B ) , (12) 

where F k is a radiometric factor related to the atmospheric transmittance 
between the target and observer and the total (direct plus diffuse) down- 
welling irradiance at the target. G k is a radiometric factor related to 
the path radiance andp^ k is the spectral reflectance of the i-th pixel for 
the k-th channel. Averaging over all pixels in a training set or for a 
model of a particular class we get 

L k (T ,u) ,p B ) = p k F k (x ,w ,p B ) + G k (x,(o,p B ). (13) 

It should be noted that the terms in Equation 13 depend upon the atmos- 
pheric parameters t and w and the background reflectance p g . If we now 
denote f(x), g(to), and h(p g ) as probability density distribution functions 
for the parameters x,u>, and p g respectively, we can then integrate each of 
these parameters over all atmospheric states for which the distribution 
functions are valid, and obtain their mean values. In a later section we 
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will perform these averages for a specific region and time of year. Thus, 
the complete average is given by 

; ■ ^ ~ F k - _G k- < i4 > 

The radiance covariance matrix is then written as 

C(L k ,L k .)«SS < L 1k ■ L k )(L ik' - 4'>' (15) 

1=1 j=l 

Substituting expressions 12 and 14 into Eguation 15 and carrying out the 
algebra gives us the final expression for the radiance covariance matrix, 

C ( Lr *Lk i ) - C(p| < ,p| < , ) C(F^,F|^i) + P |<P k 1 C ( F k , F j, i ) + 

F k F k 1 C ( p k » p k ’ ) + p k C ^ F k » G k ’ ^ + 

p k '^( F k ' »®k) + C(G k ,G k i). (16) 

Let us now examine Equation 16 carefully in order to understand its 
significance. The matrices involving F k and G k are known if we know the 
joint probability density distribution functions for the optical thick- 
nesses. For example, the covariance matrix elements for C(F k ,F k ,) are 
given by 


00 CO 


C(F k ,F k ,) -//f(F k .F kl )(F k -F k )(F k ,-F k ,)dF k dF k 


1 9 


(17) 


where f (^ k < ) is the joint probability density function. This, in turn, 
is known in terms of the distribution function g for the optical thick- 
nesses i.e., 


f ( F k 5 F k ' ) = g(V T k')/ J ’ 


(18) 
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when J is the Jacobian of the transformation between the variables. From 
reflectance models we can determine the means and covariance for the re- 
flectances in Equation 16 and, therefore, the complete covariance matrix 
for mul ti spectral radiances is determined. The diagonal value of the 
matrix is the variance. It is given by 

2 _ 2 

a 2 (L k ) = a 2 ( P k )a 2 (F 1< ) + p k o 2 (F k ) + F k a 2 (p k ) 

+ 2p k C(F k ,G k ) + a 2 (G k ) (19) 

We can now treat special cases of Equation 19. If the atmosphere is 
constant, that is, if there is no variation in the optical properties of 
the atmosphere then we have for the variance 

a 2 (L k ) = F|<a 2 (p|<)> (20) 

which indicates that as the atmospheric turbidity increases from one stable 
condition to another the variance decreases. This fact is reasonable; in- 
deed, the variance will be zero if the turbidity or the atmospheric optical 
thickness is very large because F k decreases as turbidity increases and 
there will then be no variation at all. This is illustrated in Figure 1 
for dark and bright surfaces for various atmospheres. As the visibility 
decreases the dark surfaces will brighten and the bright surfaces will 
darken and all surfaces will approach the limit point which is the path 
radiance. 

For the special case of a constant surface Equation 19 gives us 

2 

o 2 (L k ) = P k a 2 (F k ) + 2p k C(F k ,G k ) + a 2 (G k ). (21) 

Equation 21 is quite interesting; it indicates that if we know the sta- 
tistics for the atmosphere throughout a uniform spatial area then the 
reflectance can be determined by solving for P k . Thus, the atmospheric 
statistics are found by examining a data base on optical thicknesses and 
the complete variance cr 2 (L k ) is determined experimentally by remote 
sensing. A direct experiment using, say, the Landsat multispectral sensor 
should confirm the validity of Equation 21. 
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Figure 1. Two-dimensional Representation of Spectral Signatures 
as a Function of Atmospheric Turbidity 
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3.3 


Radiative Transfer 


Using an atmospheric radiation model we will now develop expressions 
for the means and variance of the path radiance. For simplicity, we will 
use a single-scattering approximation for the path radiance, i.e.. 


L p ( t ,y ,4>) 






( 22 ) 


where t r is the Rayleigh optical thickness and t is the aerosol optical 
thickness. It is assumed that the Rayleigh value for a gaseous atmosphere 
is always known. It is given by 


.(h) 


■R( h 0 ) p » 


(23) 


where x R (h Q ) is the value at some standard altitude (sealevel), P is the 
barometric pressure at altitude h, and P Q is the pressure at sealevel. The 
cosine of the nadir view angle is y and the cosine of the solar zenith 
angle is y . The factor L q is the asymptotic radiance, given by 



V 0 E o p ( y ,^,- y o ><j) o) 

47r(y+y Q ) 


(24) 


where E Q is the known extraterrestrial solar irradiance at the top of the 
atmosphere and p ( y , > _,J 0 , ) 1S 'the single-scattering phase function. The 
asymptotic value for path radiance is reached for infinite optical depth, 
i.e., for a very turbid atmosphere. A representation of the path radiance 
is given in Figure 2. We can write L 

r 

as 


L p( t) = L^(l-ae bx ) , 


(25) 
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Patti Radiance 


where 


a = e ” bt R and b = - + — . (26) 

w p o 

We now use a log-normal distribution function for the representation of the 
aerosol optical thickness i.e., 

f( T ) = — l exp — — (In T-m) 2 ] , (27) 

\/ 2 - rr ax 2a 2 J 



figure 2. Path Radiance vs Aerosol Optical Thickness 


where a 2 is the variance and m is the mean value. The mean value for Lp 
is then 


Lp = L a - L A a f e bT f(t)dT (28) 

o 


= La(I-I)- 
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Likewise, we can determine the variance in the path radiance i.e.. 


— 2 

V ( L p ) = Lp - L p , 

where 

Lp - L 2 (K-I 2 ), 

where the I and K integrals are 



and , 



(30) 


(31) 


(32) 


(33) 


where 

a - -be m and 3 o. (34) 

The integrals must be done numerically and in a later section we will eval- 
uate them to find the mean and variance in the path variance for a specific 
geographical area. 

3.4 Atmospheric Data 

The author made many contacts with atmospheric scientists and exper- 
imentalists in meteorological data analysis to acquire a comprehensive data 
base for the atmospheric optical parameters such as single-scattering 
albedo and optical thickness. Virtually nothing exists on w , so we will 
have to assume values which are generally characteristic of a particular 
region of investigation. 
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Although it is not difficult to measure the optical thickness, it 
does require spectral filters and a clock-drive to orient a sensor toward 
the sun. Usually the devices used require special handling by trained 
personnel and uniformity of devices and methods of deployment. As a re- 
sult, good values of t q (x) are available only for particular sites and for 
limited times. The most comprehensive set of turbidity values is that of 
the turbidity network or the Background Air Pollution Monitoring Network 
(BAPMON) data [5] for 1978 which can be obtained from the National Climatic 
Center in Asheville, NC. Using a two-band sensor, the data included 
turbidity measurements for selected sites on an hourly basis for the 
following locations. 


Irel and 
Spain 

Federal Republic of Germany 

Hungary 

Bulgaria 

Italy 

Turkey 

Afghanistan 

India 

Canada 

United States 
El Salvador 
Antarctica 
Pacific Islands 
Austral ia 

Data tapes are available from the National Climatic Center from which one 
can obtain daily and monthly mean values, variances, and covariances for 
the two spectral bands. 

Associated with the turbidity, one can also obtain the usual weather 
data from the National Climatic Center. Therefore, one can correlate the 
turbidity or optical thickness with various parameters such as visibility. 
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pressure, temperature, and relative humidity. Likewise, one can generate 
the multivariate probability distribution functions for all the variables. 

In addition to the above data, a gridded data base [6] is available 
for Europe containing cloud coverage, cloud type, cloud base and cloud top 
heights, as well as temperature, pressure, and relative humidity at 
eighteen layers from the surface up to 18.23 km (60,000 ft.) altitude. One 
could use this data base to generate the probability density distribution 
functions for clouds in the same way that we can do for aerosols. A multi- 
variate probability distribution function would give us the number of 
clouds within some selected area for any season as a function of cloud 
coverage, cloud type, base height, etc. Having these distributions, one 
could then calculate the expected values for the 1 i ne-of-sight transmit- 
tance for the cloud-free-1 ine-of-sight. Because we have data for various 
levels in the atmosphere, we can also generate the distribution functions 
in four-dimensions , i,e. the three space dimensions and time. Although we 
have not analyzed the cloud data, one can generate the probability distri- 
butions for an arbitrary 1 ine-of-sight for space-to-ground transmittances 
and radiances. In the present investigation we have used the turbidity 
climatology data for two channels for one month as an example of the 
technique. 

The example which we shall take is for Cerro Verde in El Salvador as 
given by the BAPMON data. This station was chosen because a fairly large 
sample of observations was' available, in comparison to some of the sample 
sizes for other months and other stations. Table 3 illustrates the aerosol 
optical thickness for the month of July, 1978 for the spectral bands 500nm 
and 880nm. The mean, variances, and standard deviations for both bands are 
calculated. Also, we determined the covariance and the correlation. It 
should be noted that the correlation coefficient of 0.979 is high which 
implies a physical connection between the two bands - certainly not a sur- 
prising result. In Figure 3 we present a histogram of the El Salvador 
data. As can be seen, most of the values lie near an optical thickness of 
0.1. In Figure 4 we have plotted the results simultaneously for both 
spectral bands. The high correlation is clearly seen. 
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Table 3 AEROSOL OPTICAL THICKNESS 


LOCATION: Cerro Verde, El Salvador 



T^feoo nm) 

1^(800 nm) 

0.19342 

0.18881 

0.18190 

0.25559 

... - -- - - - 

0.07368 

0.07368 

0.08980 

0.08059 

0.11513 

0.11513 

0.09210 

0.09901 

0.08059 

0.07599 

0.07368 

0.09441 

0.11973 

0.08980 

0.08980 

0.09441 

0.08750 

0.07368 

0.05987 

0.07599 

0.06908 

0.07829 

0.07138 

0.06447 


























































A (500 nrri) 



Covariance 


Correlation 


THICKNESS (cont.) 



x^( 800 nm) 


0.08289 


0.16348 


0.52269 


0.35920 


0.37993 


0.08059 


0.11973 






0.151 


0.01486 


0.12191 
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Similar data for the variable gaseous components such as ozone, 
water vapor and carbon dioxide are difficult to obtain. Suitable data 
bases probably exist, but to extract the effective absorber amounts and 
correlate the results with time of year and location is a large task. In 
the case of clouds, the problem is two-fold; first, one must develop 
probability distribution functions for clouds in terms of geographic 
location and time, and, second, one must then correlate the cloud type with 
effective reflectance and emissivity. A possible solution is to use 
spectral bands as in the Thematic Mapper to discriminate large, distinct 
clouds from snow. 

3.5 Results of Calculations 

Using the El Salvador data as an example, we considered the fol- 
lowing scenario: 

Solar zenith angle e = 0° 

Nadir view angle e= 0 

t r (Rayleigh optical thickness) = 0.1241 

Wavelength A= 500nm 

From these data we can determine the path radiance distribution function by 
the following: 


g ( L/Ln) = , (35) 

dL 
dr 

when f [ t ( L ) U is the probability distribution function for the optical 
thickness. The corresponding function g ( L/L A ) or g(X), where X = L/L A is 
given in Figure 5. Evaluating the integrals in Equations 32 and 33 we 
obtain 


I = 0., 40016 
K = 0.32837 
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so that the mean path radiance is 


Lp = 0. 6L^» 

and the variance in the path radiance 

V(L p ) = 0. 16824 l\. 

Similar results can be generated for 
any region and time for which optical 


(36) 
is 

(37) 

the covariance between bands and for 
thickness data are available. 



Figure 5. Probability Distribution Function for Path Radiance, 
for Cerro Verde, El Salvador (July, 1978) 
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4.0 CORRECTION ALGORITHMS 


Having determined the radiance covariance matrix and the mean 
radiance vector one can then process multispectral data for a particular 
condition. There are several methods available for processing. 


Method 1: The philosophy here is to obtain atmospheric data for large 
areas of the world for various times of day for different times of year. 
From this data base, one can then generate the statistics for the atmo- 
spheric parameters and, therefore, the corresponding probability distri- 
bution functions. Then, using radiative-transfer models one can determine 
the covariance matrices and the mean values of the radiometric parameters. 
Using actual field reflectances or models of reflectance one can determine 
the means and covariance of the surface reflectances for various classes of 
materials. Finally, all of this information is combined as in Equation 16 
to determine the actual radiance covariance matrix and the mean radiances 
for use in the decision algorithm. The resulting classification accuracy 
should be better than previous methods which rely on a "universal" mean and 
covariance. The general pattern is illustrated in Figure 6. It is obvious 
that as more atmospheric data become available the more complete will be 
our representation of atmospheric states and the more accurate will be our 
classification. 


Method 2: A variation of the above method is to consider a limited number 
of data sets. Let our reference radiometric factors be denoted by F and G, 
that is, these are averages for some training set. Let the corresponding 
averages for the unknown area be primed i.e., F 1 and G*. The bar indicates 
that we have averaged over all atmospheric states. Then, one can operate 
on the measured radiance LT to obtain a corrected radiance i.e., 



(38) 


Let us consider a simple example. From our atmospheric and material class 
data bases we determine F and G for the middle of Kansas in June. We now 
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Complete Radiance Means and 
Covariance Matrix 


Figure 6. Method for Atmospheric Correction of Multispectral Data 
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want to classsify materials in Texas in August so we use the corresponding 
values F' and G 1 for that time and location. The are then the actual 
radiances appropriate to the Texas location for each pixel i. L? is then 
the new or corrected radiance. If it turns out that the atmospheric con- 
ditions and the material classes in Kansas are the same as those in Texas, 
then IF'=F and G'=G and 

L i = 1-., , (39) 

a result which indicates that we can do no better than the training set or 
model data base. In any case, this is a method which allows us to extend 
the spectral signatures. 

Method 3: Another method we can use is based upon the assumption that one 
knows or can identify a particular class of objects in a multispectral data 
set or frame. We then have the following: 

a 2 (L) = FV(p), (40) 

where it is assumed that we know the variance a 2 ( p ) For a particular class 
of materials. From the measured radiance data we know a 2 (L). Hence we 
know the atmospheric parameter F. Likewise, in taking the average over the 
known data we have 

L = p F + G, (41) 

where L is measured and p is known. Therefore, both factors F and G are 
determined and one can process al 1 data in the region. A corrected rad- 
iance is 



where F and G are the known factors from a reference area. 


( 42 ) 
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5.0 CONCLUSIONS AND RECOMMENDATIONS 

In this investigation we have analyzed measured optical thickness 
data for one station in a simple two-band sensor network. The statistics 
were gathered for one month from the National Climatic Center in Asheville, 
North Carolina. The measurement stations are located in fourteen areas of 
the world and represent, as best as we can determine, a complete and 
comprehensive optical thickness data set. As described in this report, if 
all the data were analyzed we could generate the statistics and probability 
density distribution functions for optical thickness in terms of time of 
day, month, and geographic location. Even if this were done, it would 
represent data only for two bands in the visible part of the spectrum. 
Nevertheless, the use of the data set would be better than not using it as 
has been done in the past. 

In this work we have assumed diffuse reflectances for the targets. 
Although this is not a serious restriction, for completeness one should 
develop a corresponding formalism for non-Lambertian surfaces. This could 
have been done here, but it would have led to an unnecessarily complex 
formalism. The key development in this investigation is the expression for 
the general covariance matrix of radiances in terms of basic atmospheric 
and reflectance data. Given reflectance data and atmospheric data we can 
generate a more realistic radiance covariance matrix for particular spatial 
and temporal regions. This should lead to more accurate probabilities of 
classification than before because we are using data which more nearly 
represent the region under consideration. It should be noted that this is 
a statistical argument in that there may be occasions when one obtains a 
high classification probability because the unknown data just happen to 
possess values which are close to the reference data set. For example, if 
the reference set is Kansas in June and the unknown set is Texas in August, 
it is possible that the same atmospheric state could prevail in both places 
for the times considered. In general, however, this is unlikely, so that a 
statistical method based upon actual atmospheric data should provide the 
user with a consistently higher classification probability for a large 
number of data sets. 
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We recommend further development in this area. A complete stochas- 
tic atmospheric radiative-transfer model can now be developed i.e., one 
which includes the variation in the single-scattering albedo and the 
gaseous components. If data are lacking one could at least develop the 
mathematical formalism and use whatever theoretical models exist in atmo- 
spheric physics to establish some stochastic parameters. For actual remote 
sensing users we suggest the further analysis of optical thickness data and 
the implementation of the algorithms in this report to multispectral clas- 
sification problems. 
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